Dynamical Gaussian state transfer with quantum error correcting architecture 
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Transferring a quantum state of a light field to a memory is of particular importance. However, 
this transfer is usually hampered because the memory system is subjected to some noise and this 
can limit the performance of the state transfer to a great extent. In this paper, we consider the 
transfer of a Gaussian state of light to a linear medium memory such as an opto-mechanical oscillator 
and propose a dynamical feedback controller that suppresses the noise in the memory system. To 
protect an unknown state, the feedback scheme employs the specific configuration of the quantum 
error correction; that is, a three-mode Gaussian state having appropriate syndromes is taken as the 
input. Correspondingly, the memory consists of three independent linear systems. The syndrome 
errors are estimated continuously in time through the measurement of the output field, and the 
results are then fed back to control the system. Because the input is Gaussian and the systems are 
all linear, it is possible to formulate the problem using the framework of the celebrated classical 
Kalman filtering and linear quadratic Gaussian control. A numerical simulation demonstrates the 
effectiveness of the control scheme. 

PACS numbers: 42.50.Lc, 02.30.Yy, 42.50.Ct, 03.67.Pp 
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I. INTRODUCTION 

Transferring a quantum state of light to a memory is 
of particular importance for various purposes in quantum 
information technologies E , S SB S 0, B H [IS Ell, 
II2, 13, 14, 15, 16, 17, lj,[ial2aj2ir[2|. Candidates for a 
memory are largely divided into two categories: discrete 
variable systems such as an atom with distinct energy lev- 
els [H, i, [1 i, [5, 6, 7, 8^ i, [IS [m, [H m and continuous 
variable systems such as an opto-mechanical oscillator 
with a vibration mode [3 [H [3, [13, [H, [I^, [IS [H HI . 
Remarkably, some experimental demonstrations of quan- 
tum state transfer have been reported fll. Isl. [T6l. [isl. Il9l|. 

In reality, however, the memory performance is essen- 
tially limited by environmental noise that inevitably oc- 
curs during the transfer process. In a real experiment 
several state-of-the-art techniques should be, at least im- 
plicitly, employed to suppress such noise. However, to the 
best of our knowledge there exists no systematic method 
for suppressing the noise. Therefore, control theory for 
quantum memory needs to be explored. 

In working out this subject, a key fact is that due 
to the dynamical noise the information contained in the 
memory state is gradually erased, rather than that the 
noise suddenly vanishes it. Hence, if such an unwanted 
change of state can be monitored continuously in time, 
then the measurement result could be fed back to sup- 
press the noise. In fact, in [2], the probe light field for 
state transfer is continuously measured at the terminal of 
the optical path to indirectly detect the error during the 
transfer process. In this paper, we propose a dynamical 
control scheme that uses such continuous measurement 
results not only for the error detection but further for 
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feedback control to obtain the better performance of the 
state transfer. Fortunately, measurement-based quantum 
feedback control theory is well developed [23, 24, 25] and 
has many practical uses [H, [H [H, [^ [sO, [Sl], (32] ; hence, 
it will be also useful for our scheme. 

Here we describe two specific features of the proposed 
system-controller configuration. The first one is that the 
memory is a linear medium system such as a collective 
atomic ensemble or an opto-mechanical oscillator. In 
addition, the input state to be transferred is assumed 
to be a Gaussian state of light [s^, in which case the 
memory state also becomes Gaussian. Then, the feed- 
back control scheme has the form of celebrated linear 
quadratic Gaussian JLQG) control with Kalman filter- 
ing [2I m, IM [IqI, [In. The second feature is that the 
controller employes the schematic of quantum error cor- 
rection (QEC) for the purpose of protecting an unknown 
input state. The basic idea of QEC is to encode an un- 
known state into an enlarged codespace so that an ap- 
propriate syndrome measurement yields sufficient infor- 
mation about the error, which can then be used to cor- 
rect the error [13, [H, M, [13, [H, [H. Use of this QEC 
strategy for quantum memory means that the input is 
a three-mode unknown Gaussian state having appropri- 
ate syndromes, and it is transferred to the memory that 
is correspondingly enlarged; during the transfer process, 
the output light fields are measured continuously in time, 
yielding the estimate of the syndrome errors that can be 
fed back to control the system. 

The effectiveness of this control scheme is actually ex- 
pected from the fact that the dynamical feedback con- 
trol has successfully been applied to some QEC problems 
[IS [M [H, [H, [M [H [H. It should be noted that the 
problem considered in this paper is not a QEC problem 
itself; in QEC, the system state can be protected with 
the use of ancillary apparatus, whereas in our formalism, 
the initial state of the memory system cannot be pro- 
tected, but rather the input state to be transferred to 
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the system is what should be protected. This observa- 
tion is vital in the sense that our scheme does not violate 
the no-go theorem for Gaussian QEC [47j . 

The paper is organized as follows. In Section II, as 
a preliminary, we describe a single-mode Gaussian state 
transfer. The dynamical control scheme is addressed in 
two sections. First, the state preparation and transfer 
processes are given in Section III. Section IV explains the 
error detection scheme (Kalman filter) and the feedback 
controller (LQG controller). Finally, in Section V, we 
numerically demonstrate the efficiency of our method by 
taking, as an example of the memory, an opto-mechanical 
oscillator manipulated under ultra-low temperature. 



II. SINGLE-MODE GAUSSIAN STATE 
TRANSFER 




FIG. 1: (Color online) Opto-mechanical oscillator as an ex- 
ample of the single-mode linear quantum memory. The right 
end-mirror of the Fabry-Perot cavity serves as the oscillator, 
where a and 6 are annihilation operators corresponding to the 
intra-cavity and oscillator modes, respectively. Adiabatically 
eliminating the cavity mode a, Eq. ([T| is obtained. F denotes 
a Faraday isolator that facilitates one-way coupling between 
the input and cavity modes. 

As a preliminary, this section is devoted to describe 
a single-mode Gaussian state transfer from a light field 
to a memory, where the former is a coherent CW laser 
field and the latter is served by a linear system such as 
an opto-mechanical oscillator shown in Fig. [1] (see e.g., 
[13, [13, Ell for a detailed discussion). The input is a 
coherent state |ain)j with Sin the annihilation operator 
of this coherent light field. The reflected output mode 
Oout could be used to correct the errors occurring during 
the transfer process. In what follows we describe the 
dynamics of the memory system and show how much the 
input state jofin) is transferred to the memory. 

A general single-mode open linear dynamics is de- 
scribed by the following quantum Langevin equation: 

^ = —b - V^(ain + ao) - Vf e, (1) 

where h denotes the system annihilation operator. Here, 
V represents the coupling constant between the system 
and the input coherent light field = ain -I- oq, where 
ain is the mean and ag is the field annihilation operator. 
Moreover, we assume that the system couples with an 
unwanted single-mode noisy environment represented by 



the annihilation operator ^ with mean zero. T is the 
coupling strength. 

Let us now take the white noise approximation on the 
outer field modes oq and ^; i.e., (ao(i)aj(i')) = 5{t — t') 
and {^{t)i'^ {t')) ^ {n + l)5{t - t'), where n > repre- 
sents the strength of the noise. In the case of thermal 
noise, n is the averaged photon number. This approxi- 
mation allows us to represent the dynamics of b in terms 
of the It o-type quantum stochastic differential equation 
(QSDE) [ii,!^: 

db — i — bdt — y/va-i-ndt — ^/vdAo — vTiiS, (2) 

where A^it) = aQ{s)ds and S(i) = Jq £,{s)ds are the 
so-called quantum Wiener processes. Their infinitesimal 
increments satisfy the following quantum Ito-rule: 

dAodAl = dt, {dAof = dAl^ = dijdio = 0, 
dSdSl" = (n + l)dt, dE)dk = ndt, dE? = dE)"^ = 0. 

In this QSDE representation the input mode is written by 
dAin = a-i-ndt -\- dA^. Since the field state is the vacuum, 
we have {dAin) = aindt. Also note that (dS) = 0. 

Now, let us see the steady state of the memory system. 
Using the above Ito rule we readily obtain the ordinary 
differential equations of {bit)) and {b{t)'' b{t)) , which give 
{b{oo)) = -2y/VanJ{iy + r) and 

(Ag(oo)^) = (Ap(c^)2) = i + -i^>i, 

where q = {b + w)/^/2 and p ~ {b — ¥)/^/2i denote the 
dimensionless position and momentum operators of the 
system, respectively. This clearly illustrates that in the 
long time limit the system certainly acquires the infor- 
mation of the field input state jam) in the mean sense 
(note that ly and F are assumed to be known). However 
the fluctuation of the state could become much bigger 
than the vacuum fluctuation 1/2; in this case the state 
wrote down to the memory is far away from pure and has 
very small overlap with the input state. 

III. PREPARATION AND TRANSFER 
PROCESSES OF THE INPUT STATE 

In this paper, we study the system shown in Fig. [2] 
First, in the left box, the input state, which is a three- 
mode Gaussian state of light, is prepared. This input 
state is transferred into the memory shown in the middle 
box in the figure, that correspondingly consists of three 
identical linear systems; in the figure, opto-mechanical 
oscillators are particularly depicted as the memory sys- 
tem. Finally, the output fields are measured, and the 
results are used for feedback control. 

A notable feature of this transfer scheme is that the 
input state is allowed to contain some unknown param- 
eters. More specifically, the input state is generated by 
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FIG. 2: (Color online) A schematic of the dynamical state 
transfer with QEC architecture, where in this picture an opto- 
mechanical oscillator is taken as a memory system. 



vacuum field denoted by A2 , and one of the outputs is fur- 
tlier combined with the second ancilla squeezed vacuum 
field Art. This combination of beam splitters is called the 
iritter [39}. 

Let us describe the above encoding process in detail. 
As in the single-mode case discussed in Sec. II, we write 
the source mode as dAi^ = a-mdt + dAi, where and 
Ai are the mean amplitude and the quantum fluctuation, 
respectively. In particular, we assume 



implying that the momentum element of the mean is 
known to be zero. As mentioned above, in this case, 
the three-mode encoding is sufficient to detect the posi- 
tion errors. The quantum fluctuation Aj (j = 1, 2, 3) is a 
stochastic process satisfying (dAj) = and the following 
quantum Ito rule: 



dA,dA\ 



dAj = Mjdt, 



(Nj + l)dt, 



dAf 



dA]dAj 
M*dt, 



Njdt, 



where the parameters Nj > and Alj G C have to satisfy 
^ji-^j + 1) — l-^^iP- The fluctuation parameters Ni and 
Ml of the source mode are unknown. For the ancilla 
modes A2 and ^3, we set 



M2 = Mc, = 



embedding an unknown single-mode field into the 
three-mode light fields with the use of ancilla squeezed 
fields and beam splitters. This is no more than the 
schematic of QEC; actually, Ai^ corresponds to a source 
mode to be protected, and for this purpose it is encoded 
into a larger Hilbert space where we can construct ap- 
propriate syndromes that do not depend on the unknown 
source mode but only contain information about the er- 
rors. Those errors could be corrected by feedback con- 
trol through the syndrome measurement depicted in the 
lower box of the figure. Note that the three-mode en- 
coding allows us to detect only either the position or the 
momentum error, but the proposed scheme can be ex- 
tended to a nine-mode code to detect errors acting on 
both the position and momentum operators [s^. 

In this section, we describe the above-mentioned state 
preparation and transfer processes, and then evaluate in 
detail the memory performance, when without any error 
correction. 



A. State preparation 

For discrete variable systems, a series of CNOT gates is 
often used for encoding. A possible continuous-variable 
analogue is realized in linear optics circuit shown in the 
left box of Fig. [51 first, the information source mode Ai^ 
is mixed via a 1 : 2 beam splitter with an ancilla squeezed 



that is, the ancilla 



which satisfy Nj{Nj + 1) 

modes are the identical pure squeezed vacuum fields with 
squeezing parameter /i. Let us collect the field quadra- 
tures in a single vector as 



Wi 



(3) 



{Ql,PuQ2,P2,Q3,P3)\ 

where Qj = (i^ + i])/\/2 and Pj = {Aj - i])/V2i. 
Then, through the fritter, Wi becomes dW{ — (3dt + 
TdWi, where T is the orthogonal matrix corresponding 
to the combination of the above-mentioned 1 : 2 and half 
beam splitters: 

T = 




and 



f3 = 



V2c 



V3 



-(1,0,1,0,1,0)^. 



Using the Ito rule, it is found that the position quadra- 
tures of W{ satisfy 



d 
di 



\/i ^ j. 



(4) 
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This quantity is essentially equivalent to the power spec- 
trum density of {{Q[ — Qj)^) at the center frequency of 
the career laser field [5l|. We then find that Eq. (jU be- 
comes zero for every i,j when taking the limit ii — >■ — oo. 
This means that, in the Schrodinger picture, this ideal 
input state generated through the fritter lives in the 
codespace spanned by {\Q,Q,Q)} and is of the follow- 
ing GHZ-like form : 



m)\Q,Q,Q)dQ, 



(5) 



where ipiQ) is a Gaussian wave function of the source 
state. The tilde means that the state is the ideal one. 
Next let us focus on the following quantity: 



Vm := ((Q'l - Qi)') + ((Q'2 



3{{Pi 



Pi 



Pi 



(6) 



which leads to 

dVm 
dt 



3p/ 



Ml ~ Afi* + 1) 



If Ai,A2, and are all coherent fields, we have 
dVfd/dt = 7.5, which thus can be interpreted as the clas- 
sical limit. This means that a non-classical input state is 
generated when dVfd/dt < 7.5, which is now equivalent 
to that /i < and Mi > 0. Moreover, it is known that a 
symmetric state such as the one taken here is entangled 
if dVfd/dt < 6 [H, [5^. For instance setting a coherent 
source state, i.e., Ni = AIi ~ 0, and taking the ideal 
limit II — >■ —00, we obtain dVfd/dt — )■ 4.5; hence in this 
case the GHZ-like state ([5]) is indeed entangled. 



B. Transfer process 

Let us next describe the transfer process of the mode 
W{ to the memory served by the three identical linear 
systems. The dynamics of the memory is given by the 
combination of Eq. ([2]) with the input field replaced by 
dW{; that is, the vector of quadratures, 

with Qj — {bj + 6])/V2 and pj = (bj — Sj)/-\/2i, satisfies 
the following QSDE: 

dx = Axdt + udt - ^/T'dW[ - ^/fdW2 

= Axdt + udt - y/Tipdt + BdW, (7) 

where W = {Wi, M^a)^ and 
v + T 



A 



h, B={-^T,-Vrie 



The fluctuation vector Wi is given by Eq. ([3]) and 

W2 = V2(5R(Si),3(Si),3fJ(S2),3(S2),3?(S3),3(S3)r 



is the vector of noise processes to which the linear systems 
are subjected; hence Ej (j = 1,2,3) satisfies 



dSjdS] = {n+l)dt, dE^dEj = ndt, = dEf 



0. 



Again, n > represents the noise strength. Finally, 
u e represents the control input for the memory. Note 
here it is assumed that both the quadratures of each lin- 
ear system can be controlled. In what follows of this 
section we set u — 0. 

The state transfer can be evaluated in terms of the 
mean vector and the covariance matrix, as discussed in 
Sec. II. First, the mean vector (x) = {{qi) , . . . , (p^))^ 
obeys d{x) /dt — A{x) — ^/v/i, thus in the long-time limit 
we have (i(oo)) = V^^-^/? = -2y^l3/{v + F). That 
is, the system state becomes Gaussian with mean vector 
parallel to /?, implying that the memory certainly ac- 
quires information of the input state in the mean sense 
(note that 1/ and F are assumed to be known). Next let 
us consider the covariance matrix of the system: 

V = {AxAx'^ + {AxAx'^y)/2, Ax^x-{x). 

Using the Ito rule, the time evolution of V is found in the 
following form called the Lyapunov differential equation: 



V ^AV + VA^ + T{n + l/2)h + vTAT^ 



(8) 



where 



A = diag{Ai, A2, A3}, 
^ _ fN,+n{M,) + 1/2 5(Af,) 



5(M, 



N,-^{M,) + l/2 



Note A2 = A3 by assumption. Now, to explicitly evalu- 
ate the state transfer, let us take the source state to be 
coherent, in which case Mi — Ni = 0; then the steady 
solution of Eq. ([5]) is obtained as 



= rdiag{Fi, V2,V3}T^, = diag{z;+, «-}, 



where 

± _ v{2Nj ± 2Mj -^1) + F(l + 2n) 
~ 2{:y + r) " 

Note V2 = V3. Then we have the explicit form of the 
fidelity between the input state \ipin) and the steady state 
of the memory with its mean appropriately displaced: 



F = (l/'in|/5oo|V'in) = 

- n 



1 



(T—Q.,-\-fl,—fl 



v/dct(yoo + Vin) 

2{iy + r) 
2ve'' + T{e" -f 1 + 2n) 



(9) 



where Vm = TKT^ is the covariance matrix of the input 
state IV'in)- Note that Eq. ([5]) is the ideal limit of IV'in)- 
Eq. ^ shows that, when F = 0, we have F = \ without 
respect to the value of ^l. This means that the input 
state is perfectly transferred into the memory when it is 
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not coupled to the noisy environment. But in reality the 
system must be subjected to some noise; QEC is expected 
to overcome this issue. Motivated by some continuous 
variable QEC protocols found in the literature [s^ 113, 
Us,, .39,], let us take the operator qi — qj = 1, 2, 3) as a 
syndrome. Actually, this has desirable properties shown 
as follows. First, it can be seen from the structure of the 
system dynamics ([7]) that the syndrome does not contain 
the unknown parameters ain,A'^i, and Mi. Second, in 
the ideal limit ^ — oo, the variance of each syndrome 
is given by 



r(2n + 1) 



(10) 



which takes zero only when the system does not couple to 
the noisy environment. This means that, if in the small 
interval \t, t + dt) only one of the three linear systems 
is subjected to the noise, simultaneous "measurement" 
of the syndrome operators can detect in which system 
that error has occurred and how much it is. (The reason 
why the double quotation is taken is that the syndrome 
operator cannot be directly measured in our setting, but 
it needs to be estimated] see the next section.) Because of 
these two properties, in the limit ^ — — oo, the operators 
qi — Qj certainly plays a role of a syndrome. Here it should 
be noted that F takes the maximum value when /i = 
0; that is, without error correction the encoding process 
merely degrades the performance of the state transfer 
[s^ ]. Hence we really need the correction process. 

Before closing this section, let us evaluate the entan- 
glement of the system state, using essentially the same 
quantity as Eq. ([6]). Again in the case A^i = Mi = 0, we 
have 

Vsys ■■= (A(<?1 - q2f) + {Mh - qsf) + {A{qs - qif) 
(4.5 + 3e'')i^ (7.5 + 15n)r 



When r = and At 



-oo, it is found Vs 



4.5. 



Hence, together with (A(qi — qj) ) = in this case, the 
memory state becomes an entangled GHZ-like state. But 
in the realistic situation with F > such entanglement 
can vanish. Particularly in the ideal case fi — > — oo the 
sufhcient condition for entanglement, Vsys < 6, leads to 
n < 0.1(i^/F) — 0.1; that is, in the case of thermal noise, 
the averaged photon number of the environment field 
must be less than about one order of magnitude below 
the S/N rate v/T. 



setting, it should be appropriately estimated through in- 
direct measurement. In this section we first describe a 
continuous-time estimator, i.e., the filter, for the syn- 
dromes. We then present the dynamical feedback con- 
troller that is based on the syndrome filter. 



A. Syndrome Filter 

Let us focus on the following two kinds of operator 
vectors: 

^ f V2{pi + p2 + pa)^ 
h = Bix=^\ h+qs- 2gi | (11) 
V3{q2~q3) 



and 



P . 1 /g2 + 93 - 2^1 
VE[ V3{q2 - qs) 



(12) 



Here Bi := ZiT^ and B2 := Z2T^ are isometric matri- 
ces with 



^0 1 0> 
= ( 1 I , Z2 
,000010, 



1 
1 



If only one of the three linear systems is subjected to 
the noise, S2 can detect in which system that error has 
occurred. Hence, S2 serves as the syndrome operator. 
The first element of si is used to evaluate entanglement 
of the memory state. 

Our task is to estimate the operator si or S2 through 
certain indirect measurement. The lower box in Fig. [2] 
depicts an optical configuration that achieves this goal 
for the case of Si; that is, we construct the tritter acting 
on the output fields that are reflected at the field-system 
coupler. The observables to be measured by homodyne 
detectors with appropriate LO phase are then given by 



dYi = %/2^Jid< + ^/2{dPl, dQ2, dQ^y 



(13) 



Hence measuring Yi actually implies the indirect mea- 
surement of si. Note this contains the source fluctuation 
dPi , hence in this case only the mean value of the input 
state, /3, can be kept unknown. On the other hand, for 
the case of §2, two of the three output fields are measured; 



dY2 = V2^S2dt + V2{dQ2, dQj, 



(14) 



IV. SYNDROME FILTER AND DYNAMICAL 
FEEDBACK CONTROL 

In the previous section it was found that — qj serves 
as the syndrome through which the error can be detected. 
Since this observable cannot be measured directly in our 



Y2 does not contain both the mean and covariance of 
the source state, thus we can estimate the syndromes 
without respect to the unknown parameters. Note here 
that the position squeezing of the ancilla fields reduces 
the fluctuation dQ2 and dQa; in particular taking the 
ideal limit /i — > —00 we have di^ = ^/2vs2dt, thus this is 
no more than the syndrome measurement. 
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We now describe the filter of Si, which is constructed 
with the measurement results of (i = 1,2). Let us 
write Eq. (fT3|) or (fT4| in the following general form: 



dY = Cxdt + DdW. 



(15) 



For the case of estimating si, the matrices correspond to 
C = v^Si = v^ZiT^ and D = 72(^1,03x6)- Also 
for the case of §2, C ^ V^B2 = V^Z2T^ and D = 
\/2(^2, 02x6)- For the simple notation, we do not put the 
index i on Y, C, and D. The continuous measurement 
of Y enables us to construct the filter for the dynamics 

d7T{x) = Anixjdt + udt - ^/vfidt + Kdw, (16) 
= (V;C^ - V2^TAZ^)(2ZAZ^)-\ (17) 
dw = dy- CTr{x)dt, (18) 

where Z implies Zi or Z2. The set of above equations is 
called the quantum Kalman filter [H, HI, HI |30, iH . The 
6-dimensional c-number vector irlxt) — E(xt|3^t) repre- 
sents the quantum conditional expectation of xt condi- 
tioned on the set of measurement data = {j/s | < 
s < t} with i/g the measurement result of Yg. Note Tr{xt) 
is the least mean square estimate of xt- Vc is the condi- 
tional covariance matrix — 7t{AxAx^ + {AxAx^)^)/2 
with Ax — X — Tr{x), and it obeys the following Riccati 
differential equation: 



- 2KZAZ^K^, 



1 



(19) 



where K is given by Eq. ([T7)) . Note this is not a stochas- 
tic process. The filtering equation for the syndrome si 
or S2 can be directly obtained from Eq. as follows: 



d7r(s) = dBn{x) 

= BATT{x)dt + Budt - ^/vBI3dt + BKdw 
= A7T{s)dt + Budt + Kdw, (20) 

where A = BAB^ and K = BK. Here, B implies Bi 
or B2- Also we have dropped the index of Si and simply 
denote s. The unknown parameter ain does not appear 
in Eq. ([^0]) due to = 0. Moreover, for the case of 
estimating §2, the corresponding Kalman gain K — B2K 
does not contain Ni and il/i, thus the filter can update 
the estimate of S2 without respect to the unknown source 
state. We call Eq. (PU)) the syndrome filter. 

B. LQG control for the state transfer 

Due to the coupling to the noisy environment, the sys- 
tem state must escape from the codespace spanned by 
{k,?, but this error can be detected by estimating 
the syndrome operators. Therefore, it is expected that 



a feedback controller minimizing the estimated value of 
the syndromes corrects that error. Note that the con- 
troller must be a dynamical one so that it can deal with 
the dynamical noise; hence, the minimization should be 
carried out through the whole transfer process. Fortu- 
nately, because of the linearity of both the dynamics and 
the output equation, this requirement is satisfied when 
employing the LQG control strategy J9, 30, 31]. 

That is, we can construct an optimal control law u* that 
minimizes the following quadratic-type cost function: 



J = 



1 



u^Ru)dt^, 



(21) 



where R = > represents the penalty for the control 
input. Corresponding to the syndrome operators si or 
§2, the weighting matrix Q is set to: 

=diag{9,3,3}, Q2 =diag{3,3}. 

It is immediately seen that 

sjQiSi =(gi - 92)^ + {q2 - qzf + {q3 - Qif 
+ 3(pi +P2 +PzY, 

slQ2S2 -(91 - + ih - hf + [h - qif- (22) 

Therefore, for both cases the optimal control tries to min- 
imize the syndrome errors. In the case of si, the optimal 
control further takes into account the effect of entangle- 
ment and tries to keep the GHZ-like form of the state. 
The LQG control theory gives the explicit form of the 
(stationary) optimal controller: 



(23) 



where P is the solution to the following algebraic Riccati 
equation: 

A^P + PA + Q-PBR-^B'^P^O. (24) 

Let us here assume that R ~ tIq > 0. Then, correspond- 
ing to si and §2, the matrix P is respectively obtained 
in the following form: 



Pi =diag{/i,/2,/2}, P2=diag{/2,/2}, (25) 



where 



/! = -- 



/2 



(^ + r) 



[ly + T) 



We here give two brief remarks, (i) The first one is 
about the structure of the optimal feedback control input 
u* . In the case of si, it is given by 



3r 



/7r(g2 -qi) +7r(93 - qi)\ 
-/i7r(pi +P2 +P3)//2 
7r(gi - 92) + 7r(93 - 92) 

-/l7I'(Pl +P2 +P3)/f2 

7r(gi - 93) + T^ih - qs) 

V-/l7I"(pi +P2 +^3)7/2/ 



The first element can be regarded as a linearization to 
the following sign function (A = /2/3?'): 

A[sgn(7r(g2 - Qi)) + sgn(7r(^3 - qi))] , 

where sgn(a;) = +1 if a; > and —1 otherwise. Hence, 
in the small interval [t, t + dt), if the position shift of qi 
due to the thermal noise is larger than the others, the 
controller adds the largest inverse shift —2Xdt on dqi so 
as to cancel out the error. The other elements of u* 
have similar meanings. It should be maintained that the 
optimal controller obtained in the LQG control setup has 
this kind of natural structure found in the usual QEC 
scheme, (ii) The second remark is related to the no- 
go theorem [i^ mentioned in Sec. I. The QEC problem 
corresponding to our setting is that the system's initial 
state is to be protected by feedback control. However, 
this goal is never accomplished, because, in general, a 
stable Kalman filter forgets the initial values of 7r(io) 
and Vc(0); that is, the initial Gaussian state with mean 
7r(i;o) and covariance T4(0) cannot be protected. This 
is a way to understand the no-go theorem for Gaussian 
QEC from the dynamical control viewpoint. 
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FIG. 3: (Color online) Fidelity between the input state 
and the steady state of the controlled oscillators, versus the 
squeezing parameter and the control strength — logjr. 



V. SIMULATION 

In this section, a numerical simulation is provided to 
demonstrate effectiveness of the proposed control scheme. 
As an example of the memory we here take an opto- 
mechanical oscillator shown in FIG. [TJ The oscillator 
must be subjected to thermal noise even at ultra-low tem- 
perature; in this case the noise strength n corresponds 
to the averaged photon number n = {Qf^^-m/ksT _ 
where fc^, T, and ujm denote the Boltzmann constant, 
the temperature, and the career frequency of the ther- 
mal channel, respectively. 

The parameters are taken as follows. The oscilla- 
tor couples with the input field with strength v/2'k = 
30 kHz, while for the thermal channel it is r/27r ~ 
1 Hz. The thermal channel is with frequecy ujm/2'K = 
10 MHz and with temperature T = 4 K, which leads 
to n = 8.8 X 10'^. The mean value of the input state 
is taken as ain = —230 Hz^^^, which leads to — — 
[100,0, 100,0, 100,0]. Note that l^inp is the mean photon 
number per unit time. 



A. Coherent source state 

First let us consider the case where the source state 
is taken as a coherent state lain)- This means that the 
input fluctuation is set to A'^i = Mi = 0, implying that, 
in addition to the syndrome operators, pi -\-p2+P3 can be 
estimated. Hence, we can now apply the feedback control 
scheme based on ■k{si). 

Figure |3] shows how much the encoding and the feed- 
back control improve the fldelity between the three-mode 



input state and the oscillators' steady state. They are 
quantified as follows, (i) Encoding strength is directly 
evaluated by the squeezing parameter /i; actually, as in- 
dicated by Eq. making — /i bigger means that the 
input state is going to approximate the GHZ-like state 
([5]). (ii) The LQG optimal control input u* can take pos- 
sibly a bigger value by tuning the penalty parameter r 
smaller (Recall that we have set R — rig in Eq. ([11]).) 
Hence, — log2r is interpreted as the control strength, (iii) 
The fidelity is given by the same form as Eq. ([SJ with 
Voo now replaced by the covariance matrix of the ensem- 
ble average of the controlled system state, say V^; that 
is, F = l/^det(V4 + ^m)- Note that the controlled sys- 
tem variable Xt stochastically changes in time according 
to the QSDE ([7]) with the input u* a function of the esti- 
mate 7r(si), i.e., Eq. (^51) . Hence corresponds to the 
first 6x6 block matrix of the 9x9 covariance matrix 
of (x^, 7r(si)^)^ that satisfies the controlled Lyapunov 
equation (|A1[) ; see Appendix for the detailed calculation. 

Let us now discuss the performance of the LQG opti- 
mal feedback control. It is observed from Fig. [3] that big- 
ger control strength indicates bigger fidelity for all value 
of /J,. This means that the feedback control can always 
reduce the excess fluctuation brought from the thermal 
environment, without respect to how much the source 
state is encoded. Regarding the encoding strength, how- 
ever, care should be taken in its choice. Actually, while 
larger squeezing enables us to detect a bigger error signal 
which can induce more efficient feedback control, at the 
same time the memory state becomes more sensitive to 
the thermal noise. In other words, too much squeezing 
makes the input state fragile to the noise; consequently 
there exists an optimal value of the squeezing parame- 
ter, and it is about fj,* = —0.4. Note this value can be 
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reached within the current technology. In this case, the 
amount of improvement of the fidehty via the feedback 
control compared to the case without encoding and con- 
trol is about 0.05; but this is not a big improvement. 
One reason to this limitation is that the error taken in 
this paper is the worst one in the sense that each system 
is subjected to the thermal noise for all time; this kind of 
error is not fully tractable via the standard QEC proto- 
col. Therefore, it is expected that the proposed control 
scheme could show possibly much better performance of 
the state transfer against some weak noise such as simple 
displacement acting only one of the oscillators. 




■1 \ i 

-2 I , 

0.02 0.04 0.06 0.08 0.1 

Time [ms] 

FIG. 4: (Color online) Time evolutions of (a) 7r(g2 — ga) 
and (b) 7r(gi). For both plots, the less (red) and larger (blue) 
fluctuating lines correspond to the cases where the optimal 
feedback control is performed or not, respectively. 

Next, Fig. S] shows the time evolutions of (a) the syn- 
drome 7r((72 — 173) and (b) the first oscillator's position 
7r((7i). For both plots, the less (red) and larger (blue) 
fluctuating lines correspond to the cases where the op- 
timal feedback control is performed or not, respectively. 
The error bars show the standard deviation of the es- 
timation error, which is calculated from the solution to 
the Riccati equation ((T^ . The parameters are taken as 
^ = —0.4 and r = 10~^, which were shown in Fig. [3] 
to be the optimal values attaining the maximum fidelity. 



Note that, because a-m is unknown, qi cannot be exactly 
estimated; but we here plot the exact value of 7r(gi) just 
for demonstration. Figure S] (a) demonstrates that the 
fluctuation of the syndrome is certainly suppressed by 
the feedback control; this is the result that should be ex- 
pected, since the controller is designed to minimize the 
syndrome errors. On the other hand, it is not straightfor- 
wardly expected that controlling qi would actually work 
well, but Fig. |4] (b) demonstrates that the feedback con- 
trol based on the syndrome estimation is fairly effective 
for correcting the error acting on each position of the 
oscillator. 



B. Squeezed source state 

Here we are concerned with the situation where the 
source state is squeezed. The squeezing parameter is 
taken as fii — log(2Mi -f- 2Ni + 1) with AIi real, and 
A'^i and Mi are both unknown in addition to ain- We 
then have to use the filter for §2, which does not con- 
tain JVi and Ml . But for comparison we further consider 
a squeezed source state with known covariance and un- 
known mean value; in this case the Kalman filter for si 
can be used for feedback control. This comparison will 
reveal how much the additional information (i.e., the es- 
timate of pi +P2+P3) improves the control performance 
at the expense of limiting the class of input states. 

Figure [5] shows the fidelity between the input state 
and the controlled oscillators' steady state, versus the 
squeezing parameters /i and fii. The control penalty is 
r = 10~^. In the figure the upper and lower surfaces 
correspond to the filter for si and §2, respectively. First, 
a notable fact observed from the lower surface is that, 
when aiming to transfer the input state with the com- 
pletely unknown source state, the maximum fidelity is 
attained when fii = 0; that is, squeezing the source state 
always decreases the fidelity. This is consistent with the 
standard understanding that an unknown squeezed state 
is in general fragile to the thermal noise, because it ran- 
domly rotates the phase of the state, which as a result 
brings that state into a mixed state. On the other hand, 
the upper surface in the figure shows that, when transfer- 
ring the input state with known fluctuation, which means 
that Si can be used for estimation and control, squeezing 
the source state improves the fidelity. Note that the im- 
provement is observed if /ii > and < 0; that is, the 
source state is now momentum-squeezed while the ancilla 
states are position-squeezed. Actually in this case the in- 
put state more closely approximates the GHZ-like state, 
which is an eigenstate of pi + P2 + Ps- Therefore, the 
additional information 7r(pi + p2 + P3) should certainly 
improves the control performance of the state transfer. 
However, as in the previous case, too much squeezing 
degrades the fidelity, hence fii must be optimized. 
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FIG. 5: (Color online) Fidelity between the input state and 
the controlled oscillators' steady state, versus the squeezing 
parameters fi and /xi. The upper surface corresponds to the 
case where the estimate 7r(si) is used to control, while the 
lower surface does the case with 7r(s2) only available for esti- 
mation and control. 



VI. CONCLUSION 

In this paper we proposed a new feedback control 
scheme for state transfer, which has the form of QEC. 



In particular, due to the Gaussianity of the input state 
and the linearity of the memory dynamics, the cele- 
brated Kalman filtering and LQG feedback control were 
employed. We have considered a specific system, but 
the proposed control strategy is applicable to any linear 
quantum system, as long as the output process can be 
properly defined. Although it was shown in the exam- 
ple that the controller is not a very effective one, it is 
expected that it can show much bigger improvement for 
some linear systems having relatively weak noise. Also 
possible combination with a feedforward control scheme 
|2ll [22! could be helpful to attain better performance for 
Gaussian state transfer. 



Appendix A: Covariance matrix of the controlled 
memory state 

The whole closed-loop system of the controlled mem- 
ory and the estimator is subjected to the following QSDE: 

^{kC A- kCB^ + Bf) + {rd) 

where z — (5:^, 7r(s)^)^. Then the covariance matrix 
V; = (AzAz^-h(AzAz^)^)/2 with Az = z-{z) changes 
in time according to the following Lyapunov equation: 



dt^' " {kC a - kCB^ + bf) + i^kC A - kCB^ + bf) + (a'Z?) ^ (riT + l/2)/6 j ' 

(Al) 



Now, the covariance matrix of the controlled memory state, V ^ corresponds to the first 6x6 block matrix of Vz- In 
this case, the steady solution can be explicitly obtained as follows: 



VL=-\C^A-KC^FB)--[{A-KC^FB^ -FB) (/^) + i/2),J 

A^ ~C^K^ + B^F 



-B'F' 



{A - KC)-\FB + Ay^ + r(nT + l/2)/6 + i^TAT^j. 



r 



Note that as r — >■ (i.e., cheap control), we have — > memory state is maximally reduced down to the level of 
K (00), which is the steady solution to Eq. ([T51) . That is, estimation error, 
in this limit, the quantum fluctuation of the controlled 
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